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I.   INTRODUCTION 

THE  PROBLEM 

Air  Traffic  Control  has  become  a  safety  issue  of  great 
importance  during  the  last  decade  because  of  the  many 
near-miss  or  tragic  accidents  that  have  occurred  at  military 
and  civilian  airports  worldwide.  The  main  reasons  for  these 
accidents  are: 

*  Air  Traffic  Control  system  failure, 

*  Air  Traffic  Control  erroneous  procedures, 

*  Pilot  error, 

*  Weather  conditions, 

*  Increased  air  traffic. 

This  issue  has  even  greater  significance  at  military  air 
training  bases  because  of  the  very  high  ratio  of  student 
pilots  to  experienced  pilots  that  are  using  the  air  space 
and  because  of  the  large  volume  of  aircraft  activity  in  the 
air  in  specific  areas.  It  is  standard  military  training 
practice  for  a  large  number  of  training  aircraft  to  be 
assigned  to  the  same  radio  channel  for  aircraft-to-aircraft 
and  for  aircraft-to-air  traffic  controller  communications 
and  to  have  a  large  number  of  training  aircraft  following 
approximately  the  same  air  pattern  with  the  same  air  speed 
and  altitude.  These  conditions  can  increase  the  probability 


of  breakdowns  in  synchronization  and  communication  between 
aircraft  and  between  the  aircraft  and  the  air  traffic 
controllers  and  can  therefore  increase  the  probability  of 
accidents  occurring. 

Kalamata  Air  Force  Base  is  an  Air  Force  training  base 
located  in  southern  Greece.  This  Base  utilizes  two  types  of 
aircraft  for  training  purposes,  T-37's  and  T-2's.  The  two 
different  types  of  aircraft  used  at  Kalamata  Air  Force  Base 
ordinarily  have  different  flight  schedules,  flight 
capabilities,  and  training  missions  that  would  allow  them  to 
take-off  and  begin  their  training  flights  at  different 
times.  However  because  of  current  operational  constraints  it 
is  common  for  a  group  of  T-2  and  T-37  aircraft  to  complete 
their  missions  at  approximately  the  same  time  so  that  they 
simultaneously  return  to  the  local  Air  Traffic  Pattern.  It 
is  during  the  simultaneous  approach  of  the  returning 
aircraft  that  a  critical  safety  problem  arises  due  to  the 
traffic  congestion  and  pilot  fatigue.  Furthermore,  the 
inefficient  aircraft  schedules  interfere  with  the 
performance  of  the  scheduled  activities.  A  more  efficient 
scheduling  of  aircraft  take-offs  can  remedy  the  safety 
problems  and  also  can  diminish  the  necessity  for  aircraft  to 
wait  in  order  to  get  into  the  mission  areas. 

This  thesis  provides  a  computer  simulation  model 
programmed  in  GPSS  for  the  IBM-PC  for  analyzing  current 
procedures  for  the  efficient  management  and  control  of  the 


air  traffic  of  Kalamata  Air  Force  Base.  The  thesis  contains 
the  model  logic,  the  GPSS  program,  model  validation  and 
model  results. 


II.   MODEL  DEVELOPMENT 

A.    DESCRIPTION 

The  environment  that  is  modeled  is  the  operation  of  the 
Kalamata  Air  Force  Training  Base  in  Greece.  This  operation 
is  modeled  as  a  sequential  multiserver  limited  capacity 
queuing  system.  The  server  elements  of  the  model  consist  of 
the  base,  two  runways,  an  air  traffic  pattern  and  eight 
training  mission  areas. 

The  base,  runways  and  air  traffic  pattern  are  illustrated 
in  Figure  2.1.  The  air  traffic  pattern  is  modeled  as  a  set 
of  sequential  servers  consisting  of  two  entrance  points  (EP1 
and  EP2 ) ,  the  initial  point  (IP),  the  low  initial  point 
(LIP),  the  break  point  and  the  base  key.  EP1  is  the 
entrance  point  for  the  aircraft  returning  from  a  western 
mission  area  and  EP2  is  the  entrance  point  for  aircraft 
returning  from  an  eastern  mission  area. 

The  mission  areas  are  represented  in  the  model  as  single 
points.  Specific  mission  area  training  activity  is  not 
modeled.  The  eight  mission  areas  are  the  neighboring  areas 
around  the  airport  and  are  illustrated  in  Figure  2.2.  The 
areas  west  of  the  airport  are  numbered  one  through  four  and 
the  areas  east  of  the  airport  are  numbered  five  through 
eight. 


Figure  2.1  Air  Traffic  Pattern 


Figure  2.2  Mission  Areas 


The  calling   population   consists  of  a  basic  training 
squadron  employing  T-37  type  aircraft  an  advanced  training 
squadron  employing  T-2  type  aircraft  and  occasional 
aircraft  from  other  bases. 
B.   MODEL  LOGIC 

The  general  flow  through  the  model  for  a  typical  aircraft 
is  as  follows: 

A  take-off  time  is  scheduled  for  the  aircraft  from  one 
of  the  two  squadrons  and  the  aircraft  enters  the  base 
activity  queue  at  that  time.  If  there  is  no  runway  or  air 
space  conflict  with  other  aircraft  landing  or  taking-off, 
the  aircraft  takes-off  and  proceeds  to  the  base  departure 
point  where  it  requests  and  is  assigned  to  a  mission  area  to 
carry  out  the  scheduled  training  activities. 

If  all  training  mission  areas  are  occupied  the  aircraft 
is  assigned  to  the  mission  training  area  with  the  smallest 
waiting  line  and  in  the  case  of  ties  the  aircraft  is 
assigned  to  the  mission  area  with  the  earliest  expected 
departure  time  for  the  occupying  aircraft.  All  the  aircraft 
that  wait  for  a  mission  area  to  become  available,  maintain  a 
maximum  altitude  of  7000  feet  at  the  corresponding  areas 
and  perform  training  maneuvers  consistent  with  the  altitude 
safety  requirements.  Once  the  aircraft  enters  the  mission 
training  area  it  stays  in  the  mission  area  for  a  standard 
length  of  time  for  the  aircraft  type  and  performs  the 
scheduled  activities. 


At  the  completion  of  the  mission  area  training 
activities  the  aircraft  returns  to  the  base  and  enters  the 
returning  air  traffic  pattern  and  follows  the  sequence  of 
events  as  described  below  and  illustrated  in  Figure  2.1. 

Aircraft  returning  from  the  western  mission  areas  enter 
the  air  traffic  pattern  at  entrance  point  EP1  and  aircraft 
returning  from  the  eastern  mission  areas  enter  the  air 
traffic  pattern  at  entrance  point  EP2 .  The  sequence  of  air 
legs  in  the  base  air  traffic  landing  pattern  and  the 
aircraft  capacity  of  each  air  leg  is  as  follows. 
Air  Leg Aircraft  Capacity 

EP1-EP2  2 

EP2-IP  2 

IP-LIP  2 

LIP-Break  Point  3 

Break  Point-Base  Key  3 

Base  Key-Final  2 

If  any  of  the  first  four  air  legs  is  at  full  capacity 
when  an  aircraft  attempts  to  enter,  the  aircraft  must  orbit 
at  the  air  leg  entry  point  and  wait  for  entry.  If  the  Break 
Point-Base  Key  leg  is  at  full  capacity  when  an  aircraft 
attempts  to  enter,  the  aircraft  attempting  entry  returns  to 
the  LIP-Break  Point  air  leg  entry  point.  If  the  Base  Key- 
Final  leg  is  at  full  capacity  when  an  aircraft  attempts  to 


enter,  the  aircraft  attempting  entry  performs  "go  around", 
that  is,  it  cancels  the  landing  and  attempts  a  "close 
approach"  pattern. 

Aircraft  in  the  "close  approach"  pattern  re-enter  the  air 
traffic  pattern  at  the  Base  Key  point  if  no  other  aircraft 
is  in  the  Break  point-Base  Key  air  leg  and  if  no  other 
aircraft  is  waiting  for  take-off.  Otherwise  they  re-enter 
the  air  traffic  pattern  by  joining  the  waiting  line  at  the 
LIP  entry  point. 

Upon  being  cleared  for  landing  the  aircraft  will  either 
land  and  return  to  the  squadron  or  if  flying  constraints 
allow  the  aircraft  will  perform  a  touch  and  go  landing.  If 
a  touch  and  go  landing  is  made  the  aircraft  will  either 
enter  the  "close  approach"  pattern  or  attempt  to  re-enter 
the  air  traffic  pattern  with  associated  probabilities  of 
0.20  and  0.80.  If  there  is  a  waiting  line  at  the  LIP  point 
or  if  the  IP-LIP  air  leg  is  at  full  capacity  the  touch  and 
go  aircraft  will  either  continue  to  LIP  or  re-enter  the  air 
traffic  pattern  at  the  EP1  entry  point  with  associated 
probabilities  of  0.70  and  0.30. 

It  is  also  possible  that  the  flow  of  aircraft  in  the  air 
traffic  landing  pattern  can  be  interrupted  by  emergency 
events.  Emergency  events  such  as  engine  failure,  oil 
pressure  failure,  low  fuel,  or  landing  gear  failure,  are 
common  occurrences  that  require  the  aircraft  with  the 
emergency  to  land  as  soon  as  possible.    When  such  an 


emergency  event  occurs,  the  aircraft  with  the  emergency  pre- 
empts all  other  aircraft  in  the  landing  pattern.  During  the 
emergency  the  nonemergency  aircraft  in  the  landing  pattern 
upon  reaching  the  IP,  LIP,  or  Break  air  leg  points  leave  the 
pattern  and  return  to  the  entrance  point  EP1  and  join  the 
orbiting  waiting  line.  If  a  nonemergency  aircraft  is  on  the 
Base-Final  air  leg  when  the  emergency  situation  arises  the 
aircraft  continues  and  makes  a  full  stop  landing. 

During  the  emergency  the  aircraft  present  in  the 
mission  areas  do  not  leave  their  areas   but  after  completing 

their  scheduled  activities  orbit  at  a  lower  altitude  until 
the  emergency  ends  to  avoid  interfering  in  the  local  traffic 
pattern.  Aircraft  returning  to  the  air  traffic  pattern  will 
not  enter  the  pattern  but  will  orbit  at  the  entrance  point 
until  the  emergency  ends. 

If  an  aircraft  take  off  is  scheduled  during  an  emergency 
event  the  take-off  is  delayed  until  the  emergency  event  is 
completed. 

C.   INPUT  AND  OUTPUT 

All  input  data  for  the  model  as  well  as  suggestions 
for  model  structure  was  provided  by  the  instructor 
pilots  and  from  the  control  tower  personnel  at  the 
Kalamata  Air  Force  Base.  The  general  inputs  to  the  model 
consisted  of  individual  aircraft  characteristics  and 
performance   data,  alternate   take   off  schedules,   mission 
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area  constraints,  time  distributions  for  assignments  in  the 
mission  areas,  and  the  time  distances  between  the 
reference  points  for  the  air  legs  in  the  air  traffic 
pattern. 

The  model  outputs  contain  information  pertaining  to  the 
aircraft  of  each  sguadron,  the  air  traffic  pattern  and 
mission  area  utilization.  The  following  output  is 
available. 

1.  Aircraft  total  flight  time  distribution. 

2.  Number  of  aircraft  by  squadron  in  the  east  and 
the  west  mission  areas. 

3.  Number  of  aircraft  take-offs. 

4.  Number  of  entries  in  each  mission  area. 

5.  Number  of  entries  in  each  leg  of  the  air 
traffic  pattern. 

6.  Maximum  number  of  aircraft  waiting  in  each 
queue. 

7.  Average  wait  time  for  each  queue. 

8.  Waiting  time  distributions  for  the  following 

queues: 

a.  Air  traffic  controller. 

b.  Entry  points,  EP1  and  EP2 . 

c.  Initial  point,  IP. 

d.  Low  Initial  point,  LIP. 
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III.  GPSS  APPLICATION 

A.   BRIEF  GPSS  DESCRIPTION 

GPSS  is  the  General  Purpose  Simulation  System  language 
developed  by  IBM  for  modeling  and  simulating  queuing 
systems.  GPSS  was  used  to  model  the  Kalamata  Air  Force  Base 
air  traffic  operations.  The  GPSS  program  is  included  in 
Appendix  A. 

GPSS  uses  the  process  interaction  approach  for  modeling 
in  which  the  model  entities  are  either  temporary  or 
permanent.  The  temporary  entities  are  called 
transactions  and  the  permanent  entities  are  called 
facilities  and  storages.  The  transactions  represent  the 
calling  population  and  the  facilities  and  storages 
represent  the  service  centers.  Transactions  interact  with 
other  transactions  and  with  the  facilities  and 
storages.  In  the  Kalamata  Air  Force  Base  model  the 
calling  populations  of  aircraft  are  represented  by 
transactions  and  the  mission  areas  and  the  air  traffic 
landing  pattern  segments  are  represented  by  facilities 
and  storages. 

The  modeling  and  programming  approach  in  GPSS  is  to 
define  a  set  of  programming  statements  called  blocks 
that  represent  the  entrance  and  flow  of  the 
transactions  into  the  queuing  system    composed  of  the 
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facilities  and  storages.  There  can  be  many 
transactions  simultaneously  moving  through  the  blocks. 
At  any  point  in  time  each  transaction  is  positioned  at  a 
block  and  most  blocks  can  hold  many  transactions 
simultaneously.  The  transfer  of  a  transaction  from  one 
block  to  another  occurs  instantaneously  at  a  specific  time 
or  when  some  change  of  system  condition  occurs.  Time  in 
the  GPSS  model  is  managed  by  the  next  event  sequence 
with  the  simulation  clock  changing  at  nonuniform  discrete 
time  points  when  the  state  of  the  system  changes. 
Transactions  continue  to  move  through  the  system  until  they 
either   encounter  a  waiting  line  or  service  time  delay. 

In  GPSS  simulated  clock  time  is  an  integer  value 
whose  scale  value  is  chosen  by  the  programmer.  The  unit 
of  time  is  not  specifically  stated  but  is  implied  by 
providing  all  times  in  the  same  units.  In  the  Kalamata  Air 
Force  Base  model  the  unit  of  time   used  is  the  second. 

B.   MODEL  STRUCTURE  IN  GPSS 

A  brief  description  of  some  of  the  important 
programming  blocks  and  storage  areas  used  in  the  GPSS 
Kalamata  Air  Force  Base  model  are  contained  in  this 
section. 

GENERATE  and  TERMINATE  blocks:  Transactions  are  created 
and  enter  the  system  at  one  or  more  GENERATE  blocks  and 
are   removed   from   the   simulation  at  TERMINATE  blocks. 
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The  time  and  frequency  with  which  transactions  enter  the 
system  are  controlled  by  the  GENERATE  block.  In  the 
Kalamata  Air  Force  Base  model  GENERATE  blocks  are  used  for 
the  entry  of  training  aircraft  from  the  squadrons  for 
take  off  assignments  and  for  the  entry  of  occasional 
aircraft  from  other  bases  into  the  Kalamata  air  traffic 
landing  pattern. 

ADVANCE  block:  The  ADVANCE  block  will  hold 
transactions  for  a  specified  or  computed  number  of  time 
units.  The  purpose  of  the  ADVANCE  block  is  to  hold 
the  transactions  in  service.  In  the  Kalamata  Air  Force 
Base  model  ADVANCE  blocks  are  used  to  simulate  the  time 
delays  associated  with  take  off  delays,  training  mission 
areas  and  transit  from  point  to  point  in  the  air  leg 
segments   of   the  air  traffic  landing  pattern. 

TEST  block:  The  TEST  block  is  used  to  manage  or  transfer 
transactions  based  upon  the  test  conditions.  In  the 
Kalamata  Air  Force  Base  model  the  TEST  block  is  used  to 
prevent  aircraft  from  entering  the  system  after  the  daily 
training  period  and  to  assure  that  the  aircraft  in  the 
system  are  correctly  processed  in  order  to  complete  all 
landings  after  the  daily  training  period  ends. 

GATE  block:  The  GATE  block  is  used  as  a  gate  to 
interrupt  the  flow  of  transactions  depending  upon 
conditions  that  set  the  gate  to  "open"  or  "closed".  In  the 
Kalamata  Air  Force  Base  model  GATE   blocks   are   used   to 
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prevent  aircraft  from  continuing  in  the  air  traffic  landing 
pattern  or  the  take-off  queue  during  an  emergency   event. 

SELECT  block:  The  SELECT  block  is  used  to  direct  the 
flow  of  transactions.  In  the  Kalamata  Air  Force  Base 
model  SELECT  blocks  are  used  to  assign  the  aircraft  to  the 
training   mission   areas  after  take-off. 

JOIN,  REMOVE,  COUNT  and  MARK  blocks:  The  JOIN,  REMOVE, 
COUNT  and  MARK  blocks  are  used  to  collect,  remove, 
count  and  identify  transactions  in  the  queue.  In  the 
Kalamata  Air  Force  Base  model  if  the  training  mission  areas 
are  occupied  the  JOIN,  REMOVE,  COUNT  and  MARK  blocks 
are  used  to  determine  current  aircraft  assignments  based 
upon  the  shortest  waiting  lines  for  the  areas. 

FACILITIES  and  STORAGE  areas:  GPSS  FACILITIES  and 
STORAGE  areas  are  used  to  collect  and  hold  transactions 
for  time  delays  that  can  be  associated  with  service  or 
performance  of  the  transactions.  A  GPSS  FACILITY  can 
hold  one  transaction.  A  GPSS  STORAGE  area  can  hold  more 
then  one  transaction.  In  the  Kalamata  Air  Force  Base 
model  eight  FACILITIES  model  the  eight  mission  areas,  one 
FACILITY  models  emergency  aircraft,  one  FACILITY  models 
the  air  traffic  controller,  four  facilities  model 
aircraft  synchronization,  and  six  STORAGE  AREAS  have  been 
used  to  model  the  air  leg  segments  of  the  air  traffic 
landing  pattern  as   described   in  the  following  page. 
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FACILITY 


MODEL 


101 

102 

103 

104 

105 

106 

107 

108 

DANGER 

CNTR 

DUMY1 

DUMY2 

DUMY3 

DUMY4 


Mission  area  1. 
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C.   GATHERING  STATISTICS  WITH  GPSS 

GPSS  automatically  records  data  and  collects  queue 
statistics  for  transactions  that  pass  through  a  storage 
area.  In  addition  to  the  previously  described  GPSS 
programming  blocks  that  manage  the  flow  of  transactions 
there  are  several  block  types  that  are  specifically  designed 
to  gather  statistics  on  transactions.  These  blocks  and 
their  application  are  described  in  this  section. 

QUEUE  and  DEPART  statistic  blocks:  QUEUE  and  DEPART 
blocks  are  used  to  identify  specific  data  collection 
points  in  the  queue. 

Data  on  transactions  that  move  through  the  queue  and 
enter  and  leave  associated  QUEUE  and  DEPART  blocks  is 
collected  as  separate  sets  of  queue  statistics.  In  the 
Kalamata  Air  Force  Base  model  queue  statistics  were 
accumulated  over  the  entire  period  of  the  simulation,  in 
each  queue  of  the  air  traffic  pattern  and  in  each  queue  of 
the  mission  areas.  The  accumulated  statistics  for  each  of 
the  model  queues  are  identified  in  the  following  list. 

Queue  Statistics 


Maximum  number  of  aircraft. 

Average  number  of  aircraft 

Standard  deviation  of  the  number  of  aircraft 

Average  waiting  time. 

Number  of  entries. 
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Number  of  entries  that  did  not  wait  for  entry. 
Average  wait  time  to  enter  the  queue. 

TABULATE,  TABLE  and  QTABLE  statistic  blocks:  The 
TABULATE,  TABLE  and  QTABLE  blocks  are  used  to  collect 
data  for  frequency  and  cumulative  frequency  tables. 
The  TABLE  and  QTABLE  blocks  define  the  transaction 
characteristics  that  are  to  be  counted  and  the  range  of 
the  frequency  tables.  Data  on  characteristics  of  the 
transactions  that  pass  through  a  TABULATE  block  are 
automatically  collected.  The  TABULATE  block  may  be  used 
anywhere  in  the  GPSS  program  and  can  collect  data  on  all 
of  the  transactions  computed  characteristics.  The 
TABULATE,  TABLE  and  QTABLE  blocks  are  used  in  the  Kalamata 
Air  Force  Base  model  to  tabulate  the  total  flight  time 
distribution  for  the  two  types  of  aircraft  and  waiting  time 
distributions  for  all  of  the  model  queues. 
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IV.   MODEL  VALIDATION 

The  Kalamata  Air  Force  Base  model  was  validated  by 
comparing  model  output  to  historical  data  for  a  specific 
take  off  schedule  for  the  squadrons  of  aircraft.  Historical 
data  was  made  available  for  each  aircraft  for  a  three  day 
period.   The  historical  data  is  included  in  Appendix  C. 

For  each  aircraft  the  historical  data  consisted  of  the 
interarrival  times  for  take  off,  the  time  spent  waiting  in 
the  air  traffic  controller  queue,  the  time  spent  waiting  to 
enter  area  3  and  the  total  flight  time.  The  historical  data 
was  accumulated  and  averaged  for  the  three  days  of  base 
activity.  Using  the  historical  interarrival  rates  as  model 
inputs  the  model  generated  waiting  times  and  total  flight 
times  were  compared  to  the  historical  data. 

The  historical  interarrival  rate  appeared  to  be 
nonhomogeneous  on  a  daily  basis.  Therefore  the  operations 
day  was  broken  into  time  periods  for  which  the  interarrival 
rates  were  homogeneous.   These  time  periods  were: 

1)  07:30   -   08:40  5)  11:40   -    14:00 

2)  08:40-09:10  6)  14:00-14:30 

3)  09:10  -  11:10  7)  14:30  -  16:20 

4)  11:10  -  11:40  8)  16:20  -  17:30 
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For  each  time  period  the  interarrival  times  appeared  to 
have  an  exponential  distribution.  This  hypothesis  was 
tested  using  the  Kolmogorov-Smirnov  test  and  was  not 
rejected.  The  hypothesis  test  results  are  included  in 
appendix  D.  With  this  base  information  the  model  was  run 
for  a  period  of  eight  days  using  the  hypothesized 
exponential  distributions  for  the  take  off  interarrival 
times.   The  validation  comparison  tests  follow. 

Figure  4.1  displays  a  plot  of  the  comparison  of  the 
actual  wait  time  with  model  generated  wait  time  for  entry  to 
area  3,  and  Figure  4.1a  displays  the  regression  line  for  the 
regression  model: 

Simulated  Data=A+B*( Actual  Data). 

The  slope  of  the  regression  line  is  .97  with  a  standard 
deviation  of  .027.  The  hypothesis  B=l  was  tested  with  the 
t-test  and  was  not  rejected.  The  probability  level  for  the 
t-statistic  was  .16. 

Figure  4.2  displays  a  plot  for  the  comparison  of  actual 
wait  time  with  model  generated  wait  time  in  the  air  traffic 
controller  queue  and  Figure  4.2a  displays  the  regression 
line  for  the  regression  model: 

Simulated  Data=A+B*  (Actual  Data). 

The  slope  of  the  regression  line  is  .98  with  a  standard 
deviation  of  .043.  The  hypothesis  B=l  was  tested  with  the 
t-test  and  was  not  rejected.  The  probability  level  for  the 
t-statistic  was  .43. 
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Figure  4.3  displays  a  plot  of  the  comparison  of  actual 
total  flight  time  with  model  generated  total  flight  time 
and  Figure  4.3a  displays  the  regression  line  for  the 
regression  model: 

Simulated  Data=A+B*( Actual  Data). 

The  slope  of  the  regression  line  is  .96  with  a  standard 
deviation  of  .019.  The  hypothesis  B-l  was  tested  with  the 
t-test  and  was  not  rejected.  The  probability  level  for  the 
t-statistic  was  .09. 

These  results  provide  a  validation  of  the  model  and  show 
that  the  model  data  compares  favorably  with  the  historical 
data  for  the  specified  take  off  schedules. 
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Figure  4.1a  Regression  of  Simulated  Wait  Time 
on  Actual  Wait  Time  in  Area  3 


23 


OfiL  vs  SMftTED  UftlTING  TIMES 

AT  LINE  OP 


104 


H  34- 


p  84- 


0 


84- 


! 


54 


T 


TTT 


I    I    I    I    I    I    I    I    I   I    I    I    I    I    I    I    I    I 

I 


/■■ 

../•.' 


:/»•'■■  • 


4- 


t 


'I''' i  i  i  i  i  i  i 


0        100        200       300 


500       600 


WITMG  TIMES 


Figure  4.2   Comparison  of  Actual  Wait  Time 
with  Simulated  Wait  Time  at  Line  Up 


24 


Regression  of  SIMULATED  ON  ON,  M 
waiting  times  at  line-up  point 


'J  !    !    I    !    I    |    I    I    !    I    !    I 


95- 


II  II  I  I  I  I  I  I  I 


75 


c: 

•JJ 


v 


45  Ltiii 

44 


yy'S.y, 


_ 


* 


I  I  I 


fit        74         34         34         104 
actual 


Figure  4.2a  Regression  of  Simulated  Wait  Time 
on  Actual  Wait  Time  at  Line  Up 


25 


fiCTUfiL  vs  SIWUD  TOTAL  FLIGHT  TIME 

Aircraft  137 


100 


H  80 


60 


20 


i — rr 


y 


/ 

/• 
/• 
/• 


r 


35 


42 


i*  i'  I  »i  'i   \   i»  i1 


#  SIMULATED   - 
'■o-'flCM ~ 


45 


48 


51 


;  TIME 


J L 


54 


en 


(X 100) 


Figure  4.3  Comparison  of  Actual  Flight  Time 
with  Simulated  Flight  Time 


26 


Segression  of  Si 


TOM  FLIGHT  TIHE 


20  40 


60  80 


Figure  4.3a  Regression  of  Simulated  Flight  Time 
on  Actual  Flight  Time 


27 


V.   RESULTS-ANALYSIS 

This  section  contains  a  description  of  four  basic  and 
two  additional  take-off  schedules  for  the  squadron  of 
training  aircraft  for  Kalamata  Air  Force  Base  that  were  run 
in  the  model  and  a  comparison  analysis  of  the  resulting  air 
traffic  controller  congestion  in  order  to  determine  a 
reasonable  efficient  take-off  schedule.  In  order  to  run 
each  take-off  schedule  in  the  model,  the  GPSS  program  had  to 
be  modified.  These  program  modifications  are  included  in 
Appendix  B. 

A.  SCHEDULES 

Schedule  1:  This  schedule  consists  of  a  forty  minute 
take-off  cycle.  The  cycle  consists  of  two  consecutive 
fifteen  minute  take-off  periods  followed  by  ten  minutes  of 
no  take-off  activity.  The  T-37  aircraft  squadron  assigns 
five  aircraft  for  take-off  in  the  first  fifteen  minute 
period.  The  T-2  aircraft  squadron  assigns  five  aircraft  for 
take-off  in  the  second  fifteen  minute  period.  This  cycle  is 
repeated  until  all  the  training  aircraft  are  scheduled  for 
take-off.  Aircraft  that  can  not  take-off  in  their  assigned 
period  are  recycled  for  later  take-off.  Aircraft  taking-off 
have  priority  use  of  the  runway  over  normally  landing 
aircraft.  This  schedule  was  proposed  for  analysis  by  the 
two  squadrons  at  Kalamata  Air  Force  Base. 
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Schedule  1A:  This  schedule  is  derived  from  and  is 
identical  to  Schedule  1  except  that  the  number  of  aircraft 
scheduled  for  take-off  in  each  fifteen  minute  period  is  four 
instead  of  five.  This  schedule  was  suggested  as  a  result  of 
reviewing  the  model  output  data  for  Schedule  1.  It  was 
thought  that  this  change  would  decrease  the  number  of 
training  flights  while  also  decreasing  the  air  traffic 
congestion  and  waiting  times. 

Schedule  IB:  This  schedule  is  derived  from  and  is 
identical  to  Schedule  1  with  the  take-off  cycle  period 
extended  to  sixty  minutes.  This  schedule  was  also  suggested 
as  a  result  of  reviewing  the  model  output  data  for  Schedule 
1.  It  was  thought  that  this  change  would  also  decrease  the 
number  of  training  flights  while  decreasing  the  air  traffic 
congestion  and  waiting  times. 

Schedule  2:  This  schedule  consists  of  an  eight-minute 
take-off  cycle.  The  cycle  consists  of  two  consecutive 
four-minute  take-off  periods.  The  T-37  aircraft  squadron 
assigns  one  aircraft  for  take-off  in  the  first  four-minute 
period.  The  T-2  aircraft  squadron  assigns  one  aircraft  for 
take-off  in  the  second  four-minute  period.  This  cycle  is 
repeated  until  all  the  training  aircraft  are  scheduled  for 
take-off.  Aircraft  that  can  not  take-off  in  their  assigned 
four-minute  period  are  recycled  for  later  take-off.  This 
schedule  was  also  suggested  as  a  result  of  reviewing  the 
model  output  data  for  Schedule  1.   It  was  thought  that  this 
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model  output  data  for  Schedule  1.  It  was  thought  that  this 
change  would  increase  the  number  of  training  flights  while 
decreasing  the  air  traffic  congestion  and  waiting  times. 

Schedule  2A:  This  schedule  is  derived  from  Schedule  2  and 
is  identical  to  Schedule  2  except  that  the  take-off  cycle  is 
extended  to  ten  minutes  with  two  five-minute  take-  off 
periods.  This  schedule  was  suggested  as  a  result  of 
reviewing  the  model  output  data  for  Schedule  2. 

Schedule  3:  This  schedule  is  a  reproduction  of  the 
schedule  currently  in  use  at  Kalamata  Air  Force  Base.  The 
schedule  contains  no  structure  and  take-off  times  are 
scheduled  at  random.  The  distributions  of  current 
interarrival  times  for  take-off  times  were  analyzed  and 
found  to  be  exponential  for  different  periods  during  the 
day.  This  analysis  is  contained  in  Appendix  C  and  was  also 
used  for  model  validation. 

B.    RESULTS 

The  model  output  data  for  each  schedule  are  contained  in 
Appendix  E.  The  measures  of  effectiveness  used  in  comparing 
schedules  were  performance  and  efficiency/safety. 
Performance  is  measured  by  the  average  number  of  training 
aircraft  scheduled.  Efficiency  is  measured  by  the  average 
waiting  times  in  the  mission  area  queues.  Safety  is 
directly  related  to  efficiency  in  that  the  smaller  waiting 
times  mean  less  hazardous  flying  conditions.   The   empirical 
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distributions  for  both  of  these  measures  for  each  schedule 
are  contained  in  Appendix  E. 

The  summary  results  for  Schedule  1  are  of  interest  because 
this  schedule  was  suggested  by  the  training  squadrons  at 
Kalamata  Air  Force  Base.  This  schedule,  by  comparison  with 
the  other  schedule  results,  does  not  provide  the  highest 
values  for  performance  and  efficiency/safety.  It  was 
because  of  these  results  that  the  other  schedules  were 
derived  from  Schedule  1  by  sensitivity  analysis. 

The  summary  results  for  Schedule  3  are  of  interest 
because  this  schedule  is  a  reproduction  of  the  schedule 
currently  in  use  at  Kalamata  Air  Force  Base.  This  schedule 
by  comparison  with  the  other  schedule  results,  also  does  not 
provide  the  best  values  for  performance  and 
efficiency/ safety. 

It  was  found  that  by  applying  Schedule  1A  and 
comparing  the  results  to  Schedule  1  that  performance 
decreased  but  that  efficiency/safety  improved  drastically. 
Figure  5.1  compares  the  cumulative  mission  area  waiting  time 
distributions  and  Figure  5.2  compares  the  cumulative  LIP 
point  waiting  time  distributions  for  Schedules  1  and  1A. 

Schedule  2  results  were  more  efficient  than  Schedule  1 
results.   However,  by  applying  Schedule  2A  mission  area 
waiting  times  decreased  dramatically  and  provided   the  most 
efficiency/safety  measure  for  all  schedules  examined. 
Figure  5.3  compares  the  cumulative  mission  area  waiting  time 

31 


distributions  and  figure  5.4  compares  the  cumulative  LIP 
point  waiting  time  distributions  for  schedules  2  and  2A. 
Figure  5.5  compares  the  cumulative  mission  area  waiting  time 
distributions  for  schedules  1A,  2A  and  3.  It  is  obvious 
that  Schedule  2A  is  preferable  to  Schedule  1. 

Schedule  3  results  were  more  efficient  than  Schedule  1 
results.  The  performance  measure  for  Schedule  3  however  is 
the  lowest  of  all  the  schedules  examined.  Figure  5.6 
compares  the  cumulative  LIP  point  waiting  time  distributions 
for  schedules  1A,  IB,  2A  and  3. 

C.   ANALYSIS 

From  the  above  summary  results  the  preferred  schedules 
appear  to  be  Schedules  1A  and  2A.  An  analysis  of  variance 
was  performed,  using  the  function  "ANOVA"  from  the  OA3660 
APL  WORKSPACE,  to  test  the  hypothesis  that  the  mean  mission 
area  waiting  time  differences  for  schedules  1,  1A,  2,  2A  are 
not  significant.  The  analysis  of  variance  results  are 
contained  in  Table  5.1.  These  results  show  that  the  null 
hypothesis  of  no  significant  differences  between  mean 
mission  area  waiting  times  is  rejected  at  significance 
levels  greater  than  .995. 

The  Sum  of  Squares  from  the  previous  analysis  was  broken 
into  three  components  in  order  to  test  for  individual 
effects  rather  than  just  a  schedule  effect  using  individual 
degrees  of  freedom.    The  results  of  this  analysis  are 
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contained  in  the  tables  5.2  and  5.2a.  These  results  show 
that  the  hypothesis  of  no  significant  difference  between  the 
compared  mean  mission  area  waiting  times  for  the  selected 
schedules  is  rejected  for  each  comparison  at  significance 
levels  greater  than  .975. 

TABLE  5.1 
ANALYSIS  OF  VARIANCE  RESULTS 


ANOVA  TABLE 

SOURCE       DF 

SS          MS        F 

SCHEDULE      3 

232864.74    77621.59    27.77 

ERROR         4 

11181.05     2795.27 

TOTAL         7 

244045.78 

R-SQUARE  =  0.954 

OVERALL  MEAN  = 

200.68 

TREATMENT  EFFECTS 

-87. 56266. 68***. 50  14.39 

-68.5     24.02 

0.435  -17.92 

68.5    -24.02 

-0.435   17.92 
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TABLE  5.2 
ANOVA  WITH  INDIVIDUAL  DEGREES  OF  FREEDOM 
SCHEDULE  1       SCHEDULE  1A     SCHEDULE  2A    SCHEDULE  2 


44.62 

181.62 

49.138 

443.34 

7.62 

6.75 

197.15 

232.99 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

-1 

-1 

-1 

-1 

1 

1 

-1 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

-1 

1 

1 

1 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-1 

D  =  SK ■  •  2 

b  l  j 

Z=SbijXj 

Wi^Z^/D. 
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64162.575 

4 

-708.48 
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4 

415.77 
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2 

-137.00 

9384.500 

2 

48.04 

1153.920 

2 

0.87 

0.378 

2 

1284.50 
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TABLE  5.2a 
ANALYSIS  OF  VARIANCE  RESULTS 
USING  INDIVIDUAL  DEGREES  OF  FREEDOM 


SCHEDULES 

SS 

DF 

MS 

F 

(l,lA)vs(2,2A) 

1  vs  1A 

2  vs  2A 

RESIDUAL 

64162.575 

125485.970 

43216.173 

11181.048 

1 

1 
1 
4 

64162.575 

125485.970 

43216.173 

2795.262 

22.954 
44.892 
15.460 

TOTAL 

244045.760 

7 
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Figure  5.4  Cumulative  LIP  Waiting  Times 
Schedule  2  and  2A 
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VI .   CONCLUSIONS  -  SUMMARY 

A.   CONCLUSIONS 

If  an  added  emphasis  is  to  be  placed  upon  the 
efficiency/safety  factor  in  scheduling  the  aircraft  then 
Schedule  2A  is  the  preferred  schedule.  If  Schedule  2A 
conflicts  with  local  base  operations  due  to  other  operating 
constraints  then  Schedule  1A  is  the  next  best  schedule.  The 
difference  in  mean  mission  area  waiting  times  for  these  two 
schedules  is  made  more  significant  when  it  is  also 
realized  that  landing  aircraft  have  runway  priority  in 
schedule  2A  and  do  not  in  Schedule  1A.  Figure  6.1  compares 
the  cumulative  mission  area  waiting  time  distributions  and 
figure  6.2  compares  the  cumulative  LIP  point  waiting  time 
distributions  for  schedules  1A  and  2A. 

If  added  emphasis  is  to  be  placed  upon  the  performance 
factor  in  scheduling  the  aircraft  then  Schedule  2  is  the 
preferred  schedule.  However,  the  increase  of  twenty  four 
scheduled  aircraft  is  at  the  expense  of  more  than  a  ten  fold 
increase  in  mean  mission  area  waiting  times  that  contributes 
to  air  traffic  congestion  and  pilot  fatigue  and  should  be 
avoided  if  possible.  Again,  if  Schedule  2  conflicts  with 
local  base  operations  due  to  other  operating  constraints 
then  Schedule   1   is  the   next   best   schedule.    It  is 
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recommended  that  the  efficiency/safety  factor  be  the 
deciding  factor  in  selecting  schedules. 

Table  6.1  contains  a  summary  comparison  between  the 
different  schedules.  In  this  table  the  first  column 
contains  the  average  number  of  flights,  the  second  column 
contains  the  percentage  of  the  aircraft  that  waited  in  the 
mission  area  queues,  the  third  column  contains  the 
percentage  of  the  aircraft  that  waited  more  than  180  seconds 
in  the  mission  area  queues,  the  fourth  column  contains  the 
percentage  of  the  aircraft  that  waited  more  than  30  seconds 
at  LIP  and  the  fifth  column  contains  the  average  conditional 
waiting  time,  in  seconds,  in  the  mission  area  queues. 

TABLE  6.1 
SUMMARY  COMPARISON  BETWEEN  SCHEDULES 


#  of  A/C 
scheduled 

%  of  A/c 
waiting 

in 
m.  areas 

%  of  A/C 
waiting 
>180"in 
m. areas 

%  of  A/C 
waiting 
>30"in 
LIP 

Avg . Conditional 
waiting  time 
in  seconds 

Sch.2A 

106 

2 

0 

0.8 

11 

Sch.lA 

106 

7 

5.6 

1.6 

45 

Sch.2 

130 

26 

18 

2.5 

203 

Sch.lB 

89 

7 

7 

2 

230 

Sch.3 

99 

14 

13 

3 

368 

Sch.l 

130 

26 

17 

3.5 

496 
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Figure  6.3  compares  the  cumulative  mission  area  waiting 
time  distributions  and  figure  6.4  compares  the  cumulative 
LIP  point  waiting  time  distributions  for  schedules  1  and  2. 
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B .   SUMMARY 

This  thesis  presents  an  IBM-PC  GPSS  model  for  studying 
aircraft  take-off  schedules  at  Kalamata  Air  Force  Base  in 
Greece.  The  model  is  specific  to  Kalamata  Air  Force  Base 
because  of  the  structure  and  sequence  of  the  queues  used  in 
modeling  the  air  traffic  control  points  and  the  training 
mission  areas.  The  GPSS  program  can  be  transferred  to  any 
IBM-PC  compatible  computer  that  will  run  GPSS.  The  model 
can  therefore  be  used  at  the  Kalamata  Air  Force  Base  to 
continue  to  examine  other  aircraft  take-off  schedules  and 
to  help  reducing  fuel  consumption. 

Six  aircraft  take-off  schedules  were  examined  and  a 
comparison  of  results  was  based  upon  factors  of  performance 
and  efficiency/safety.  The  take-off  schedule  currently  used 
at  Kalamata  Air  Force  Base  is  shown  to  be  a  poor  performer 
with  a  low  efficiency/safety  factor.  The  take-off  schedule 
proposed  for  use  by  the  training  squadron  personnel  at 
Kalamata  Air  Force  Base  is  shown  to  have  the  highest 
performance  factor  and  the  worst  efficiency/safety  factor. 

The  specific  schedule  of  ten  minutes  of  a  ten  minute  take- 
off cycle  consisting  of  two  consecutive  five  minute  take-off 
periods  from  each  squadron  with  runway  priority  given  to 
landing  aircraft  is  shown  to  be  the  best  schedule  based  upon 
the  recommended  emphasis  placed  upon  efficiency/safety 
rather  than  on  numerical  performance. 
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APPENDIX   A 

This    appendix    contains    the    GPSS    program    for    the    Kaiamata 
Air    Force    Base    model. 
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; CHECK   FLIGHT   TIME   FOR  T-2  A/C 
;G0   FOR   FULL   STOP 

;20%   GO   TO    IPE,    OTHERS   CONTINUE 

;807.   GO   TO    IPE 

■IF   NO   A/C   FOR   TAKE   OFF   CONTINUE 

; DOING   CLOSE   PATTERN 

;BASE   KEY 

; CHECK   THE   QUEUES    IN   ALL   AREAS 

'; CHECK    IF   ALL   QUEUES   ARE   OCCUPIED 

;FIND   THE   EARLIEST  OCCUPIED   AREA 
,AREA      ; REMOVE   FROM   THE   GROUP 

; SELECT  AREA  WITH   THE   MIN  QUEUE 

,AREA      ; REMOVE   FROM   THE   GROUP 

;A/C   NOT   GOING   TO   MISSION   AREAS 
•PERFORM   SHEDULED   ACTIVITIES 
;BACK   TO  THE  AIRPORT 


CAP1 

CAP2 

CAP3 

CAP4 

AREA1 

AREA2 

AREA3 


QTABLE 
QTABLE 
QTABLE 
QTABLE 

QTABLE 
QTABLE 

QTABLE 


580   AREA4      QTABLE 


MAKE        TABLES 

MP7, 3600, 240, 18 

MP6, 3600, 240, 18 

CNTR, 0,60, 40 

CAP1,0,30,40 

CAP2,0,30,40 

CAP3,0,30,40 

CAP4,0,30,40 

101,0,120,32 
102,0,120,32 

103,0,120,32 
104,0,120,32 


; TOTAL  TIME  DISTRIBUTION 
TOTAL  TIME  DISTRIBUTION 
TIME  DISTRIBUTION  IN  CNTR  QUEUE 
TIME  DISTRIBUTION  IN  CAP1  QUEUE 
TIME  DISTRIBUTION  IN  CAP2  QUEUE 
TIME  DISTRIBUTION  IN  CAP3  QUEUE 
TIME  DISTRIBUTION  IN  CAP4  QUEUE 
;TIME  DISTRIBUTION    IN  AREA  Q101 

;TIME   DISTRIBUTION      IN   AREA  C102 
;TIME  DISTRIBUTION    IN   AREA   Q103 

;TIME   DISTRIBUTION      IN   AREA   Q104 


2  "■  9  -" 
2400 


2430 

2440 
2450 
2460 

2470 
2480 
2490 
2500 

2520 


S40 


AREAS 
AREA6 
AREA? 

AREAS 

FIRST 

* 

* 

CAP1 
CAP2 
CAP3 
CAP4 
CAPS 
CAP6 
* 

* 

BNRT 


BVARIABLE 

DEFINE 

STORAGE 
STORAGE 
STORAGE 
STORAGE 

STORAGE 
STORAGE 


10/  ,0,  ^.- ,■-._ 
1  ("iP  ■"  1  ^>m  t*: 


I  line  Uib  i  KibwJ  i  iUN 
;TIME  DISTRIBUTION 

,   1  .  i _  UlDlnlDUl  iUn 


A  l"s 

AREA 

IN 

AREA 

IN 

AREA 

IN 

AnlEA 

(P2 ' E ' 101+P2 ' E ' 102+P2 ' E ' 103+P2 ' E ' 104 ) 
THE        STORAGES 


2 
2 

7 


TIME   SEGMENT 


iJDv 

2570 


GENERATE 
TEST  E 
TERMINATE 


N*D0WN,N$FINISH 

1 


;TIME  ARRIVES  AT  30600 
;WAIT  UNTIL  ALL  THE  A/C  LAND 
j SHUT  OFF  THE  RUN 
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This    appendix    contains    the    GPSS    program    modification, 
for    the    different    take    off    schedules. 

SCHEDULE    1 


10   ***********»**#****#***#*************************^ 

20  *  * 

30   *  AIR-TRAFFIC   CONTROL      SIMULATION  * 

40   *  * 

50   *  * 

60   *******************************************^ 
70   * 

90  *  EXPONENTIAL  DISTRIBUTION 

ion  * 

110  EXP     FUNCTION     RN1.C24 

0,0/ . 1 , . 1 04/ . 2 , . 222/ . 3 , . 355/ . 4 , . 509/ . 5 , . 69/ . 6 , . 915/ . 7 . 1 . 2/ . 75 , 1 . 38 


." §,l!  6/1  84,  l!  83/7e3,2~  12/79, 2.  3/.  92, 2.52/. 94,2. 81/.  95. 2. 99/. 96, 3. 2 

.97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995. 5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 

120  * 

130  * 

140  MEAN    FUNCTION     C1,D8 

3600 . 380 . 1 1 /5400 , 22 1 . 60/ 1 3200 , 483 . 62 

1 5000 , 242 . 63/23400 , 557 . 83/25200 , 124 . 58 

3 1 800 , 494 . 60/36600 , 22 1 . 30 

150  * 

160  * 

170         GENERATE 

180  SWITCH!  LOGIC  S 

190         ADVANCE 

200         LOGIC  R 

210         ADVANCE 

220         TRANSFER 

230         GENERATE 

240  SWITCH2  LOGIC  S 


250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 

460 

470 
480 


ADVANCE 

LOGIC  R 

ADVANCE 

TRANSFER 

GENERATE 

SEIZE 

ADVANCE 

RELEASE 

TERMINATE 

GENERATE 
ASSIGN 
MARK 
TRANSFER 


l"1 

300 

2100 
,SWITCH1 
900  1 

300 

2100 

.SWITCH2 

9000,FN$EXP 


CREATE  A  SINGLE  TRAN3AC" 
SET  LOGIC  SWITCH  1 
T-37  TAKE  OFF 
RESET  LOGIC  SWICH  1 
•NONE  T-37  TAKES  OFF 


600,180 
DANGE 


9100,FN$EXP 


4 
,EPI 


;SET  LOGIC  SWICH  2 
:T-2  TAKE  OFF 
; RESET  LOGIC  SWICH  2 
;NONE  T-2  TAKES  OFF 

; CREATE  EMERGENCY  EVENTS 

; EMERGENCY  HOLDS 

; EMERGENCY  TERMINATES 

:A/C  FROM  OTHER  BASES 
•.ASSIGN  TO  A  PARAMETER 
;MARK  THE  TIME 


CREATE  AIRCRAFTS  FROM  TWO  DIFFERENT  SQUADROMS  ON  BASE 


GENERATE 
TEST  L 
GATE  LS 
ASSIGN 

TRANSFER 

GENERATE 


60,10 

N$GNRT,1,PISTA 

1,PISTA 

1.1 

,  DOWN 

60,10 


; CREATE  T-37'S 

;THE  SIMULATION  STOPS  AFTER 

; ASSIGN  TO  A  PARAMETER 


: CREATE  T-2'S 
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SCHEDULE  li 


10  it**************************************************************** ***++ 
20  *  * 

*  AIR-TRAFFIC  CONTROL   SIMULATION  * 

*  * 

*  + 

* 
* 

*  EXPONENTIAL  DISTRIBUTION 


RN1.C24 

~i/ .  4 , .  509/ .  5 , .  69/ .  6 , .  915/ .  7 , 1 .  2/ .  75 , 


30 

40 

50 

60 

70 

80 

90 

ion 

110  EXP     FUNCTION 

0 , 0  / .  1 , .  1 04  / .  2 , .  222  / .._.,.  -jjj/ 

. 8 , 1 . 6/ . 84 , 1 . 83/ . 88 , 2 . 12/ . 9,2.3/ . 92 ,2. 52/ . 94 ,2. 81 / .95,2. 99/ .96,3. 

. 97 , 3 . 5 / . 98 , 3 . 9 / . 99 , 4 . 6/ . 995 , 5 . 3/ . 998 , 6 . 2/ . 999 , 7/ . 9998 , 8 

120  * 

17.0  * 

140  MEAN    FUNCTION     C1.D8 

3600 , 380 . 1 1 /5400 , 22 1 . 60/ 1 3200 , 488 . 62 

1 5000 , 242 . 63/23400 , 557 . 83/25200 , 1 24 . 5S 

31 800 , 494 . 60/36600 , 22 


: CREATE  A  SINGLE  TRANSACTION 
;SET  LOGIC  SWITCH  1 
;T-37  TAKE  O^F 
•RESET  LOGIC  BWICH  1 
;NONE  T-37  TAKES  OFF 


150  * 

160  * 

170 

GENERATE 

» » » i 

180  SWITCH1 

LOGIC  S 

i  ■  ' 

190 

ADVANCE 

7.  or. 

200 

LOGIC  R 

1 

210 

ADVANCE 

2100 

220 

TRANSFER 

, SWITCH! 

230 

GENERATE 

V 900,1 

240  BWITCH2 

LOGIC  S 

250 

ADVANCE 

300 

260 

LOGIC  R 

9 

270 

ADVANCE 

2100 

280 

TRANSFER 

.SWITCH2 
9000,FN*EXP 

290 

GENERATE 

300 

SEIZE 

DANGE 

310 

ADVANCE 

600.180 
DANGE 

TOT) 

RELEASE 

370 

TERMINATE 

340  * 

350 

GENERATE 

9100,FN*EXP 

360 

ASSIGN 

1,2 

370 

MARK 

4 

380 

TRANSFER 

,EPI 

390  ♦ 

400  *   CREATE  AIRCRAFTS 

FROM  TWO  DIFFEI 

410  * 

420 

GENERATE 

75,10 

N$GNRT,1,PISTA 

1,PISTA 

430 

TEST  L 

440 

GATE  LS 

450 

ASSIGN 

1,1 

460 

TRANSFER 

,DOWN 

470  * 

480 

GENERATE 

75,10 
N$GNRT,1,PISTA 

490 

TEST  L 

500 

GATE  LS 

2,PISTA 

510 

ASSIGN 

10 

520  DOWN 

QUEUE 

CNTR 

530 

GATE  FV 

DANSE 

540  * 

GATE  SE 

CAPfc 

;SET  LOGIC  SWICH  2 
:T-2  TAKE  OFF 
; RESET  LOGIC  SWICH  2 
;NONE  T-2  TAKES  OFF 

; CREATE  EMERGENCY  EVENTS 

; EMERGENCY  HOLDS 

; EMERGENCY  TERMINATES 

;A/C  FROM  OTHER  BASES 
; ASSIGN  TO  A  PARAMETER 
•MARK  THE  TIME 


; CREATE  T-37'S 

;THE  SIMULATION  STOPS  AFTER 

; ASSIGN  TO  A  PARAMETER 


; CREATE  T-2'S 


; ASSIGN  TO  A  PARAMETER 
;WAIT  OUT  OF  THE  RUNWAY 

•DONT  MOVE  IF  EMERGENCY  HOLDS 
; CHECK  FOR  A\C  ON  BASE  LEG 


SCHhOiJLE 


+  ■*■ 

*  AIR-TRAFFIC   CONTROL      SIMULATION  * 

*  * 

*  * 
* 


EXPONENTIAL  DISTRIBUTION 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100  * 

110  EXP     FUNCTION     RN1,C24 

0 , 0/ .  1 , .  104/ .  2 , .  222/ . 3 , . 355/ . 4 , . 509/ . 5 ,. 69/ . 6, . 915/ . 7 , 1 . 2/ . 75 , 1 . 3 

.8, 1.6/. 84, 1.83/. 88, 2. 12/ . 9, 2. 3/. 92, 2. 52/. 94, 2. 81/. 95, 2. 99/. 96, 3. 

.97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995, 5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 

120  * 

130  * 

140  MEAN    FUNCTION     C1,D8 

3600 , 380 .  1 1 /5400 , 22 1 . 60 / 1 3200 , 438 . 62 

1 5000 , 242 . 63/23400 , 557 . 83/25200 , 1 24 . 58 

31800,494.60/36600,221  ' 


CREATE  A  SINGLE  TRANSACTION 
SET  LOGIC  SWITCH  1 
T-37  TAKE  OFF 
RESET  LOGIC  SWICH  1 
NONE  T-37  TAKES  OFF 


;SET  LOGIC  SWICH  2 
;T-2  TAKE  OFF 
•RESET  LOGIC  SWICH  2 
;NONE  T-2  TAKES  OFF 


150  * 

160  * 

170 

GENERATE 

180  SWITCH1 

LOGIC  S 

190 

ADVANCE 

300 

200 

LOGIC  R 

i 

210 

ADVANCE 

3300 

220 

TRANSFER 

,SWITCH1 

230 

GENERATE 

v 900,1 

240  SWITCH2 

LOGIC  S 

250 

ADVANCE 

300 

260 

LOGIC  R 

2 

270 

ADVANCE 

37.00 

280 

TRANSFER 

.SWITCH2 

9000,FN*EXP 

DANGE 

290 

GENERATE 

300 

SEIZE 

310 

ADVANCE 

600,180 
DANGE 

320 

RELEASE 

330 

TERMINATE 

7.40  * 

350 

GENERATE 

9100,FN$EXP 

360 

ASSIGN 

1,2 

370 

MARK 

4 

380 

TRANSFER 

,EPI 

7.90  * 

400  *   CREATE  AIRCRAFTS 

FROM  TWO  DIFFEI 

410  * 

420 

GENERATE 

60,10 

N$GNRT,1,PISTA 
1,PISTA 

430 

TEST  L 

440 

GATE  LS 

450 

ASSIGN 

1,1 

460 

TRANSFER 

,DOWN 

470  * 

480 

GENERATE 

60.10 
N$GNRT,1,PISTA 

490 

TEST  L 

500 

GATE  LS 

2,PISTA 

510 

ASSIGN 

1,0 

520  DOWN 

QUEUE 

CNTR 

530 

GATE  FV 

DANGE 

540  * 

GATE  SE 

CAP6 

; CREATE  EMERGENCY  EVENTS 

; EMERGENCY  HOLDS 

; EMERGENCY  TERMINATES 

;A/C  FROM  OTHER  BASES 
■ASSIGN  TO  A  PARAMETER 
;MARK  THE  TIME 


; CREATE  T-37'S 

;THE  SIMULATION  STOPS  AFTER 

; ASSIGN  TO  A  PARAMETER 


; CREATE  T-2'S 


ASSIGN  TO  A  PARAMETER 
WAIT  OUT  OF  THE  RUNWAY 
DONT  MOVE  IF  EMERGENCY  HOLDS 
CHECK  FOR  A\C  ON  BASE  LEE 
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SCHEDULE   2 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100   * 

110   EXP 

0,0/. 1, 


********************************************************************** 

*  * 

*  AIR-TRAFFIC   CONTROL      SIMULATION  * 

*  * 

*  * 
********************************************************************** 
* 
* 
* 


EXPONENTIAL   DISTRIBUTION 


FUNCTION 
104/. 2, .222/. 3 


RN1,C24 
'55/. 4,. 509/. 5,. 69/. 6, 


V|  Vf    ■    *  ■  ■  i  \F~TJ    m  4tm  m  a  *-+—£- i    »•_'■•  'JJU/    ■  T  a  a  JV   f  f    a   vJ  ■    iU//    iUi   i7l  J/    i    .■     %    1   •   Z.  /   >    /  «_J  ,    a   .   -J'O 

.8, 1.6/. 84, 1.83/. 88, 2. 12/ . 9. 2. 3/. 92. 2. 52/. 94, 2. 81/. 95, 2. 99/. 96, 3. 2 

.97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995, 5. 3/. 998, 6. 2/. 999, 11. 9998, 8 

120   * 

17.0    * 

140   MEAN  FUNCTION  C1,D8 


3600 , 380 . 1 1 /540C , 22 1 . 60 / 1 3200 , 488 . 62 
15000 , 242 . 63/23400 , 557. 83/25200 , 124 .  58 
3 1 800 , 494 . 60 / 36600 , 22 1 .  30 
150   * 


160   * 

290 

GENERATE 

9000, FN* EXP 

300 

SEIZE 

DANGER 

310 

ADVANCE 

600,180 

320 

RELEASE 

DANGER 

330 

TERMINATE 

340   * 

350 

GENERATE 

9100,FN*EXP 

360 

ASSIGN 

370 

MARK 

4 

380 

TRANSFER 

,T0EP1 

790   * 

400   *        CREATE   AIRCRAFTS 

FROM   TWO   DIFFEI 

410   * 

411 

GENERATE 

120, 10, ,4 

412 

ASSIGN 

.DOWN 
120,10,60,4 

413 

TRANSFER 

414 

GENERATE 

415 

ASSIGN 

,D0WN 

416 

TRANSFER 

417   * 

420 

GENERATE 

480,60,1500 

430 

TEST   L 

N*GNRT,1,PISTA 

431    * 

450 

ASSIGN 

,D0WN 

460 

TRANSFER 

470    * 

480 

GENERATE 

480,60,1740 

490 

TEST   L 

N*GNRT,1,PISTA 

491    * 

510 

ASSIGN 

1,0 

520   DOWN 

QUEUE 

CNTR 

530 

GATE  FV 

DANGER 

540 

GATE   SE 

CAP6 

550 

SEIZE 

CNTR 

560 

DEPART 

CNTR 

561    * 

562  * 

L    I 

N  E        UP 

570   LINEUP 

ADVANCE 

70,20 
CNTR 

580 

RELEASE 

590 

MARK 

6 

600 

MARK 

7 

CREATE   EMERGENCY   EVENTS 


EMERGENCY   TERMINATES 

A/C  FROM  OTHER  BASES 
ASSIGN  TO   A  PARAMETER 
MARK   THE   TIME 


SOUADROMS   ON   BASE 

; CREATE   THE   FIRST   4   T-37'S 
; ASSIGN  TO  A  PARAMETER 

; CREATE   THE   FIRST   4  T-2'S 
; ASSIGN   TO  A   PARAMETER 


; CREATE   T-37'S 
;THE   SIMULATION   STOPS   AFTER 
:ALL   THE   A/C   HAVE   LANDED 
; ASSIGN   TO   A   PARAMETER 


CREATE  T-2'S 

THE  SIMULATION   STOPS   AFTER- 
ALL  THE  A/C  HAVE   LANDED 
ASSIGN  TO  A  PARAMETER 
WAIT  OUT  OF  THE  RUNWAY 
DONT  MOVE    IF   EMERGENCY   HOLDS 
CHECK  FOR  A\C  ON   BASE  LEG 
CAPTURE   THE   CONTROLLER 
GOING  FOR  LINE   UP 


LINE-UP   CHECK 
TAKE-OFF 

START   FLIGHT   TIMEFOP   T-2   A/C 
START   FLIGHT   TIME   FOP    T-3"7   A/C 
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SCHEDULE  2A 


10  ***+**********♦*************************■*-***************************** 
T)  *  * 

30  *  AIR-TRAFFIC  CONTROL  SIMULATION  * 

40  *  * 

50  *  * 

60  *************************************************************** 
70  * 
80  -& 

90  *  EXPONENTIAL  DISTRIBUTION 

100  * 

110  EXF  FUNCTION     RN1,C24 

0, 0/ . 1 , . 104/ . 2 , . 222/ . 3 , . 355/ . 4 , . 509/ . 5 ,. 69/ . 6 . . 915/ . 7 , 1 . 2/ . 75 , 1 . 38 
.8, 1.6/. 84, 1.83/. 88, 2. 12/. 9, 2. 3/. 92, 2. 52/, 94, 2. 81/. 95, 2. 99/. 96, 3. 2 
.  97 , 3 . 5/ . 98 , 3 . 9 / . 99 , 4 . 6 / . 995 , 5 . 3/ . 99S , 6 . 2/ . 999 , 7/ . 9998 , 8 
120  * 
130  * 

140  MEAN    FUNCTION     C1.D8 
3600 ,  330 . 1 1  /5400 ,  22 1 .  60/ 1 3200 ,  488 .  62 
15000 , 242. 63/23400 ,557. 83/25200 , 124. 58 
31800 , 494. 60/36600 , 221 . 30 
150  * 
* 

; CREATE  EMERGENCY  EVENTS 

; EMERGENCY  HOLDS 

; EMERGENCY  TERMINATES 

;A/C  FROM  OTHER  BASES 
; ASSIGN  TO  A  PARAMETER 
;MARK  THE  TIME 


160  * 

290 

GENERATE 

9000,-N*EXF 

300 

SEIZE 

DANGER 

310 

ADVANCE 

600,180 

DANGER 

-T<-!/-v 

RELEASE 

330 

TERMINATE 

340  * 

350 

GENERATE 

9100,FN$EXP 

360 

ASSIGN 

1,2 

370 

MARK 

4 

380 

TRANSFER 

,T0EP1 

390  * 

400  *    CREATE  AIRCRAFTS 

FROM  TWO  DIFFEI 

410  * 

411 

GENERATE 

120, 10,, 4 

412 

ASSIGN 

11 

,66wn 

120,10,60,4 

413 

TRANSFER 

414 

GENERATE 

415 

ASSIGN 

,  DOWN 

416 

TRANSFER 

417  * 

420 

GENERATE 

600,60,1500 

430 

TEST  L 

N$GNRT,1,PISTA 

431  * 

450 

ASSIGN 

,5own 

460 

TRANSFER 

470  * 

480 

GENERATE 

600,60,1800 

490 

TEST  L 

N*GNRT,1,PISTA 

491  * 

510 

ASSIGN 

1,0 

520  DOWN 

QUEUE 

CNTR 

530 

GATE  FV 

DANGER 

540 

GATE  SE 

CAP6 

550 

SEIZE 

CNTR 

560 

DEPART 

CNTR 

561  * 

562  * 

L  I 

N  E   UP 

570  LINEUP 

ADVANCE 

70,20 

580 

RELEASE 

CNTR 

590 

MARK 

6 

600 

MARK 

7 

; CREATE  THE  FIRST  4  T-37'S 
; ASSIGN  TO  A  PARAMETER 

; CREATE  THE  FIRST  4  T-2'S 
; ASSIGN  TO  A  PARAMETER 


•CREATE  T-37'S 
;THE  SIMULATION  STOPS  AFTER 
;ALL  THE  A/C  HAVE  LANDED 
•ASSIGN  TO  A  PARAMETER 


CREATE  T— 2 ' S 

THE  SIMULATION  STOPS  AFTER 
ALL  THE  A/C  HAVE  LANDED 
ASSIGN  TO  A  PARAMETER 
WAIT  OUT  OF  THE  RUNWAY 
DONT  MOVE  IF  EMERGENCY  HOLDS 
CHECK  FOR  A\C  ON  BASE  LEG 
CAPTURE  THE  CONTROLLER 
GOING  FOR  LINE  UP 


LINE-UP  CHECK 
TAKE-OFF 

START  FLIGHT 
START  FLIGHT 


riMEFOR  T-2  A/C 
riME  FOR  T-37  A/C 
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SCHEDULE   3 


20   *  * 

30   *                                        AIR-TRAFFIC   CONTROL      SIMULATION  * 

40   *  * 

50   *  * 

60  ****************************#**#*»*#**#**#****#»****#»***###*#***##^#* 

70  * 

80  * 

90   ♦  EXPONENTIAL   DISTRIBUTION 

100  * 

110   EXP  FUNCTION  RN1,C24 

0,0/. 1,. 104/. 2,. 222/. 3,. 355/. 4,. 509/. 5,. 69/. 6.. 915/. 7, 1.2/. 75, 1.38 

. 8 ,1 . 6/ . 84 , 1 . 83/ . 88 , 2. 12/ . 9,2.3/ . 92 ,2. 52/ . 94 ,2. 81 / . 95,2. 99/ . 96 , 3. 2 

.97, 3. 5/. 98, 3. 9/. 99, 4. 6/. 995, 5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 

120  * 

130  * 

140   MEAN 


FUNCTION 


C1.D8 

7nA 


3600 , 380 . 1 1/5400 , 22 1 . 60/ 1 3200 , 488 . 62 

15000 , 242. 63/23400 , 557 . 83/252UC , 124. 58 

31800,494. 60/36600,221 . 30 

150  * 

160   * 

290 

300 

310 

320 

330 

340  * 

350 

360 

370 

380 

390   * 

400    * 

410  * 

420 

430 

435 

450 

460 

470   * 


BENERATE 

SEIZE 

ADVANCE 

RELEASE 

TERMINATE 

GENERATE 
ASSIGN 
MARK 
TRANSFER 


9000.FNSEXF 
DANGER 
600,180 
DANGER 


9100,FN*EXP 

I'2 

,T0EP1 


; CREATE   EMERGENCY   EVENTS 

; EMERGENCY   HOLDS 

; EMERGENCY   TERMINATES 

;A/C   FROM   OTHER   BASES 
; ASSIGN   TO  A  PARAMETER 
;MARK  THE   TIME 


CREATE   AIRCRAFTS   FROM   TWC   DIFFERENT   SQUADROMS   ON   BASE 


510   ACFTT2 

520  DOWN 

530 

540 

550 

560 

570   LINUP 

580 

590 

600 

610 

615 

620   * 

630 

640   AREA 

650  * 

660  * 

670  * 

680 

690 

700   WAI" 


GENERATE  FN$MEAN,FN$EXP 

TEST   L  N$GNRT,I,PISTA 

TRANSFER  .49,,ACnT2 

ASSIGN  1,1 

TRANSFER  ,DOWN 

ASSIGN  1,0 

QUEUE  CNTR 

GATE   FV  DANGE 

GATE   SE  CAP6 

SEIZE  CNTR 

DEPART  CNTR 

ADVANCE  70,20 

RELEASE  CNTR 

MARK  6 

MARK  7 

ADVANCE  100,10 

TRANSFER  .25,, OTHER 

SELECT   E  2,101,108,0,F,QUEUP 

TEST   E  BV$FIRST,I,WEST 

EAST        AREAS 


SAVEVALUE 
SAVEVALUE 
QUEUE 


13+, PI 
14+, 1 
P2 


; CREATE   AIRCRAFT 

•THE   SIMULATION   STOPS   AFTER 

: ASSIGN   TO  A  PARAMETER 


ASSIGN   TO   A   PARAMETER 

WAIT   OUT   OF   THE   RUNWAY 

DONT   MOVE    IF   EMERGENCY   HOLDS 

CHECK   FOR   A\C   ON   BASE   LEG 

CAPTURE   THE  CONTROLLER 

GOING  FOR   LINE   UP 

LINE-UP   CHECK 

TAKE-OFF 

START  FLIGHT   TIMEFOR   T-2  A/C 

START   FLIGHT   TIME  FOR   T-37   A/C 

AFTER   TAKE   OFF  TO  DEPART   POINT 

FORMATIONS,    INSTR.FL. ,    CPM 

:FIND   EMPTY   AREA    IF  EXISTS   ONE 
; CHOOSE  EAST  OR  WEST   AREA 


; RECORDS   T-37'S 

; RECORDS   TOTAL   #   OF   A/C 

;WAIT    IN   THE   QUEUE 
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APPENDIX    C 

This   appendix  contains   the  historical   data   for   a  three  day 
period.       For   each   aircraft   the   historical    data   consisted  of   the 
interarrival    times    for    take    off,     the    time    spent    waiting    in    th< 
air    traffic     controller    queue,     the    time    spent    waiting    to    entei 
area   3    and   the   total    flight   time   of   T-37    aircraft.      All   times   an 
stated   in  seconds. 


Interarrival   times   for  take  off 


INARRDAY1     {Interarrival  Tme5  for  take  off  Day  1) 
557  213  254  399  968  46  16  91  211  490  484  622  43  11  7  541  41  171  160  177  191  435 
472  62  27  646  305  1410  154  490  196  479  358  1637  207  84  568  548  1211  511 
31  225  65  565  41  149  776  279  157  362  575  1147  305  288  475  363  635  118  644 
274  1623  117  173  6  31  58  24  57  183  775  92  7  142  260  7  461  77  260  587  77 
39  376  225  235  766  70  163  200  79  322  22  59  118  298  70  3  280  146 


INARRDAY2    (Interarrival  Tites  for  take  off  Day  2) 
291  228  259  597  89  70  351  1415  147  179  1006  388  167  43  478  22  135  18  21  357  257 
606  207  413  17  260  95s?  414  182  1424  486  33  213  304  1168  362  261  177  211 
91  53  310  14  128  267  304  90  37  274  11  804  959  580  425  1491  280  887  6  235 
947  114  155  141  157  46  63  152  155  151  113  126  37  1101  1142  919  242  132 
1460  216  397  6i  1C29  393  26  221  80  52  534  85  99  261  382  60  112 


INARRDAY3    (Interarrival  Tises  for  take  off  Day  3) 
43B  170  1353  577  124  144  e8  338  260  636  330  27  547  15  427  52  127  42  267  23  11 
"4  323  141  596  115  796  421  148  171  657  214  20  371  106  1356  78  120  22  51 
114  64  '22  1135  481  1832  665  296  45  34  519  230  204  14  102  708  1238  1134 
1027  570  564  916  221  49  333  7  168  87  91  2  504  252  35  62  39  55  144  165  67 
107  2  116  146  357  270  980  1422  869  334  1097  873  277  447  312  621  274  30  21 
492  71  452  530  741  143  676 
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Time  spent  waiting  in  the  air  traffic  controller  queue 


Waiting 
Time 


0 
60 
120 
180 
240 
300 
420 
480 
540 
660 
720 
840 
960 
1020 


Control 1 

er  queue 

dayl  day2  day3 

42 

45 

49 

42 

38 

31 

12 

20 

11 

3 

5 

4 

2 

2 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

Total  flight  time  aircraft  T-37 

Time  dayl  day2  day3 

3960 
4080 
4200 
4320 
4440 
4560 
4800 


6 

5 

8 

23 

18 

13 

12 

23 

17 

10 

9 
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4 

3 

2 

0 

2 

1 

0 

1 

0 

Area  3 

dayl  day2  day3 

14 

12    13 
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0     1 

1 

1     0 

0 

1     0 

0 

0     1 

0 

1     0 

0 

0     0 

0 

0     1 

0 

0     0 

1 

0     0 

1 

0     0 

0 

0     1 

1 

0     0 

0 

1     0 
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APPENDIX  D 

This  appendix  contains  the  summary  statistics  and  the 
hypothesis  test  results  for  the  eight  time  periods  for  which 
the  interarrival  times  for  take  off  were  homogeneous.  These 
times  periods  were: 

1)  07:30  -  08:40  5)   11:40  -  14:00 

2)  08:40-09:10  6)   14:00-14:30 

3)  09:10  -  11:10  7)   14:30  -  16:20 

4)  11:10  -  11:40  8)   16:20  -  17:30 

The  Tables  D.l  and  D.2  contain  the  summary  statistics 
for  each  time  period  and  the  figures  D.l  through  D.8  contain 
the  distribution  fitting  and  the  Kolmogorov-Smirnov  test 
results. 
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Table  D.l  Summary  Statistics 


Variable 


INTERVALl 


INTERVAL2 


Sample  si^e 

Aver age 

Med  ian 

Mode 

G e ometric  m e a r ■ 

Variance 

S t a  r  i  d  a r d  d  evi a  t i o  r  i 

Standard  error 

Min  imam 

Maximum 

Range 

L_  o  w  e  r  qua  r  t  i  1  e 

Upper  quartile 

I  n  i:  e  r  q  u  a  r  t  i  1  e  r  a  n  g  e 

Skewnesa 

S t a ndardized  skewnes s 

Kurtosis 

S t  a  ndardiz e d  k u r  t  o s i  s 


35 

45 

380 .114 

221 .689 

260 

160 

259 

11 

248.624 

116.221 

120438 

44135.7 

347.041 

210. 085 

58 . 6607 

31.3176 

16 

7 

1415 

798 

1399 

791 

144 

42 

547 

413 

403 

371 

1.65105 

0.928257 

3 . 98766 

2.54214 

2.56701 

-0.118849 

3.09996 

-0. 162741 

Var  .tab  1  e 


INTERVALS 


INTERVAL4 


Cample  si;:e 

Average 

Median 

Mode 

Geometric  mean 

Variance 

S  t a n  d  a  r d  d  e  v  i  a  t  i  on 

G  t  an  d  a  r d  e  r  r  o  r 

MJ nimum 

Ma: imum 

Range 

L  o w  e  r  qua  r  t  i 1 e 

U  p  pe  r  quartile 

Interquartile  range 

SI-:  ewntiss 

Standardized  skewness 

Kurtosis 

Standardized  kurtosis 


47 

488. 

617 

305 

304 

284. 

774 

229164 

478. 

711 

69. 

8271 

20 

1832 

1812 

120 

646 

526 

1 

.28888 

3 

60734 

0 

.696499 

0 

974683 

27 
242.63 
204 

14 

144 

49612 


42 

14 

778 

764 

45 

304 

259 

1 

0 

o 


641 


738 
866 


1451 
42911 
434659 
461025 
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Table  D.2  Summary  Statistics 


Var  i ab 1 e : 


INTERVALS 


INTERVAL 6 


Sample?  si 

z  a 

36 

39 

Average 

557 .833 

124.897 

Median 

4  50 

107 

Made? 

383 

*"5 

Geometric 

:  mean 

322 . 289 

69.1587 

Variance 

190050 

20045.8 

Standard 

devia t 

ion 

435-947 

141.583 

Standard 

error 

72.6579 

22.6714 

Minimum 

6 

*"> 

Ma;-;  lmum 

1623 

775 

R  a  n  g  e 

1617 

773 

Lower  que 

u~ 1 1 .1  e 

228 

39 

Upper  qui 

trtile 

901.5 

155 

Interq u. a r t lie  r 

ange 

673 .  5 

116 

Skewness 

0.705837 

3. 12258 

S  t  a  n  d  a  r  d  i  z  e  d  s  k 

ewnes 

s 

1.72894 

7.96105 

Kurtosis 

-0. 301508 

12.1983 

Standard :. 

.zed  kurtosi 

s 

-0 . 36927 

15.5498 

Variable: 


INTERVAL 7 


INTERVALS 


Sample  size 

Average 

Median 

Mode 

G  e  a  m  e  t  r  i  c  m  e  a  n 

Variance 

S  t  a  n  d  a  r  d  d  e  v  i  a  1 1  o  n 

Standard  error 

M j  n i mum 

Ma;<  imum 

Range 

Lower  quartile 

U p per  q u a r t i 1 e 

Interquartile  range 

Skewness 

S  tandard i  zed  skewness 

Kurtosis 

Standardized  kurtosis 


38 
494.605 
345 .  5 

77 
327.136 
165600 

406 . 94 

66.0143 

26 

1460 

1434 

216 

869 

653 

0 


■0 


90161 

269 
289484 


0.364259 


30 

130.5 
112 

129.0.11 
42845 
206 . 99 
37.7911 

741 

738 
70 

1.10583 
2.4727 
0.246429 
0.275516 
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Frequency  Histogram 
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Estimated  KOLMOGOROV  statistic  DPLU5  =  0.0733633 
Estimated  KOLMOGOROV  statistic  DMINUS  =  0.0923791 
Estimated  overall  statistic  DN  =  0.0923791 
Approximate  significance  level  =  1 


Figure  D.l   Distribution  Fitting  Interval  1 
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Estimated  KOLMOGQROV  statistic  DPLUS  =  0.134797 
Estimated  KGLMOGOROV  statistic  DMINUS  =  0.111456 
Estimated  overall  statistic  DN  =  0.134797 
Approximate  significance  level  =  0.386891 


Figure  D.2   Distribution  Fitting  Interval  2 
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Estimated  overall  statistic  DN  =  0.0722221 
Approximate  significance  level  =  1 


Figure  D.3   Distribution  Fitting  Interval  3 
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Estimated  KOLMOGOROV  statistic  DPLUB  =  0.0935019 
Estimated  KOLMOGOROV  statistic  DMINUB  =  0.0871474 
Estimated  overall  statistic  DN  =  0.0935019 
Approximate  significance  level  =  1 


Figure  D.4   Distribution  Fitting  Interval  4 
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Frequency  Histogram 
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Figure  D.5   Distribution  Fitting  Interval  5 
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Estimated  KOLMOGOROV  statistic  DPLUS  =  0.128466 
Estimated  KOLMOGOROV  statistic  DMINUS  =  0.0882622 

Estimated  overall  statistic  DN  =  0.128466 
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Figure  D.6   Distribution  Fitting  Interval  6 
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Figure  D.7   Distribution  Fitting  Interval  7 
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Estimated  KGLMOGOROV  statistic  DPLUS  =  0.0867169 

Estimated  KOLMOGOROV  statistic  DMINUG  ==  0.0760762 

Estimated  overall  statistic  DN  =  0.0867169 

Approximate  significance  level  =  1 


Figure  D.8   Distribution  Fitting  Interval  8 
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APPENDIX    E 

This    appendix    contains    the    GPSS    prograir.    outputs     lor 
each    take    off    schedule. 

SCHEDULE     1 


FACILITY 


ENTRIES   UTIL. 


AVE.  TIME  AVAILABLE   OWNER  PEND  INTER  RETRY  DELAY 
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0 
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"Ar'LE 

CNTR 
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STD.DEV.  RETRY  RANSE 
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TABLE  tfEAN      SIT.DEV,  RETRY  RANGE 

AREA;  ^ 1 .  92 


AREA7 
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STD.DEV. 

RETRY 

RANGE 
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1A 
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105 
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13 
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